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1. INTRODUCTION 



In the same work as he proposed the existence of a large reservoir of comets in the outskirts of the solar system, 
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ABSTRACT 

We examine the effects of passing field stars on the angular momentum of a nearly radial orbit of an 
Oort cloud comet bound to the Sun. We derive the probability density function (PDF) of the change in 
angular momentum from one stellar encounter, assuming a uniform and isotropic field of perturbers. 
' We show that the total angular momentum follows a Levy flight, and determine its distribution 

, function. If there is an asymmetry in the directional distribution of perturber velocities, the marginal 

probability distribution of each component of the angular momentum vector can be different. The 
O ' constant torque attributed to Galactic tides arises from a non-cancellation of perturbations with an 

impact parameter of order the semimajor axis of the comet. When the close encounters are rare, 
the angular momentum is best modeled by the stochastic growth of stellar encounters. If trajectories 
passing between the comet and sun occur frequently, the angular momentum exhibits the coherent 
growth attributed to the Galactic tides. 

Qh ' Subject headings: comets: general — Oort Cloud — solar system: formation 

2 ■ lOortl ((1950) suggested a two-stage process for the creation of such a cloud. First, perturbations from the planets 
increase the semimajor axes of nearby smaller objects. These interactions leave the periapses of the small bodies in 
^ ' the planetary region, but will eventually deliver enough energy to eject them from the solar system. The second stage 
of Oort cloud formation requires that perturbations external to the solar system deliver angular momentum to the 
comets. This raises their periapses out of the realm of planetary influence and saves them from eventual ejection. 
Further perturbations are necessary to lower their periapses again so that they can return to the planetary region and 
be observed from Earth. 

Oort's original suggestion for both the circularization and delivery mechanisms is the influence of other stars in 
^ _ the Galaxy as they encounter the solar system. Each star that passes the solar system delivers a small kick to each 
. comet that depends on the mass of the star, its velocity, and its distance of closest approach. iHeisler fc Tremaind 
• ' ([1986) explored the effects of a large scale planar symmetry in the swarm of stellar perturbers to find a smooth 
\ torque similar in magnitude to or even dominant over the stochastic stellar perturbations. This effect is known as 
the "Galactic tidal torque," since it can be attributed to the gradient of the average potential of the Galactic disk. 
Several groups have used numerical sim ulations to investigate the formation of the Oort cloud from a combination of 
stellar perturbations and Galactic tides (jDuncan et al.lll987l : iDones et al.ll2004f ): in all cases the two effects have been 
implemented separately. 

Recent studies have provided analytic solutions to several other stocha s tic scattering p r oblem s that arise in orbital 
dynamics and planet formation ([Collins fc Saril[2006t ICollins et al.ll2007[ ). ICollins fc Saril (|2008[ ) investigated the evo- 
lution of an initially circular orbit interacting impulsively with unbound perturbers. They showed that the probability 
■ - - per unit time of perturbing a circular orbit to an eccentricity of order e is proportional to e~^. This power law is enough 
to determine that t he eccentricity of the binary diffuses as a Levy process, and the scale of the distribution grows 
linearly with time ([Shlesinger et allll995| ). Such evolution is fast compared to the common Brownian motion- type 
diffusion where the distribution evolves as only the square-root of time. 

In t his work, we apply the framework developed for perturbations around nearly circular orbits in ICollins fc Saril 
([lOOSl), hereafter CS08, to the case of perturbations around nearly radial orbits. Section [2] presents the effects of a 
single stellar passage on a zero angular momentum comet. In Section |3| we derive and solve the Boltzmann equation 
that describes the accumulation of the changes in angular momentum from an isotropic distribution of perturbers. 
Section |4| describes the perturbations that arise from an anisotropic velocity distribution, and explains the connection 
between the effects of stellar encounters and the tidal force from the Galactic potential. Section [5] summarizes our 
conclusions. 

2. A SINGLE STELLAR PASSAGE 

In this section we discuss the change in angular momentum, eccentricity, and periapse of a comet on a nearly radial 
orbit. We call the central body of the system the "sun." We denote the position of the comet as rf,(<), and its velocity 
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vi,{t). We write the magnitude of rf,(i) as ri,{t) — \ri,{t)\, and the unit vector as Hit) = rb{t)/ri,{t). Since a radial 
orbit is by definition a straight fine, fbit) is constant in time. Furthermore, the direction of the velocity, Vbit), is 
either aligned or anti-aligned with ff,. The orbital energy per unit mass of the comet, £, sets the semimajor axis, a, 
and the orbital period, Torb- The angular momentum vector, J is zero, and the eccentricity vector is then given by 
e — Wh X 3 / (GMq) — fb — — fb- Finally, determining the position of the comet as a function of time requires specifying 
the time that the comet passes through periapse, r. 

We call each perturber a "star," and write the velocity of the star Vp. The mass of the star, rup, will typically be 
about the same magnitude as Mq, the mass of the sun; both are very large compared to the mass of the comet, rric- We 
focus this analysis on the regime where the path of the star is unaffected by the gravity of the Sun, or GMq / (bvp) <^ 1. 
Then the position of the perturber as a function of time is given by rp{t) = h + Vp{t — to) , where b describes the 
closest position of the star relative to the sun and to is the time at which the star reaches this position. 

We consider, at first, encounters between the star and the sun that occur with b ^ 2a, such that the perturbation 
to the sun-comet system can be treated in the tidal limit. We will show that these interactions are important for 
setting the angular momentum distribution when it is near zero. In section 4 we derive the evolution of the angular 
momentum as it evolves under all types of encounters including b < a. 

In this case the tidal acceleration as a function of time is given by: 



arit') = Gnipnit') 



Vpih ■ Vp) - h (h + Vpt){h ■ fb) 



(62 + {vpt')^)y^ (62 + {Vpt')^)5/^ 



(1) 



where we have translated the time coordinate by to to simplify the expression. 

The acceleration caused by each passing star affects the shape of the comet's orbit. Since the angular momentum 
is initially zero, the small impulses have a large relative effect on J. In contrast, single perturbations to e and a are 
always small compared to their initial magnitudes. The periapse of the comet, which is important for determining 
the influence of the planets on the comet, is related to the angular momentum, J — ^j2GMQq. For these reasons, we 
focus on understanding the effects of the stellar perturbations on the angular momentum vector. 

To find the total change in angular momentum for one stellar passage, we integrate the acceleration over the motion 
of the star and of the comet: AJ = J Yb{t') x aT{t')dt' . There are two limiting cases where we can evaluate this integral 
to find a closed form solution. The first is the impulsive regime, where b/vp <^ Torb- The comet spends most of its 
time with rb ~ a; however, for rare interactions that occur when rb{t) <^ a, impulsiveness requires b/vp <^ rb{t)/vb{t). 
We treat the comet as stationary over the duration of an impulsive perturbation: rb{t') — rb{to), and find the change 
in angular momentum to be: 



AJ(io) 



2Gmprb{to)'^ 
Vpb"^ 



ih X Vp){fb ■ Vp) ~ 2{fb X b){fb ■ b) 



(2) 



The other simplifying case is a very non- impulsive encounter {b/vp ^ Torb). When each orbit is very short relative 
to the timescale of the perturbation, the acceleration at each point along the perturber's path is experienced by the 
entire span of the comet's orbit. The disparate timescales in this regime allows the integral over the motion of the 
comet to be separated from the integral over the path of the star. The result is a AJ that is independent of to and r: 



AJ = ^^IHp ^ Vp){fb ■ Vp) - 2{fb X b){fb ■ b) 



(3) 



This result also follows from replacing rb{to)'^ in Equation [2] with its time averaged value, (r^) — (5/2)a2 



3. LEVY FLIGHT BEHAVIOR 

Successive perturbations cause the angular momentum delivered to the comet to accumulate. Individual pertur- 
bations add to the existing angular momentum vectorially: Jncw = x (vf, + Au) = J + AJ. Holding fb constant 
restricts the angular momentum vector to a plane. We accordingly treat J as a two-dimensional vector throughout 
this work. Since the perturbations by passing stars occur randomly, we employ the same statistical approach as CS08. 
We study the evolution of J by deriving a distribution function, f{J,t), that specifies the probability that the comet 
will have an angular momentum within the region (fj around J at time t. If we assume that the perturbations occur 
isotropically, there is no preferred direction for the accumulated angular momentum of the comet. We then expect 
that f{J,t) = f{J,t). The probability of finding the comet's angular momentum with a magnitude between J and dJ 
in any direction is 2iTf{J,t)JdJ. We relax the assumptions of isotropy in Section|4l 

We express the probability density function (PDF) for single perturbations as a frequency per unit angular mo- 
mentum, TZ{ J'). This function describes the probability per unit time that the comet receives a perturbation with 
a magnitude between J' and J' + dJ'. Given the properties of the ensemble of perturbing stars, we compute the 
frequency with the following expression: 



Tl{J') J S{\A3{vp,h,to,mp)\ - J')J^{vp,mp)vpS{h ■ Vp)d^hd^Vpdmpd{to /To 



b • 



(4) 
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where the function J^(vp, rrip) is the combined phase space density of perturbers in Vp and nip, normahzed such that the 
total mass density of perturbers in real space is p = J mpj-'{vp, mp)d^Vpdmp. This equation is analogous to Equation 
9 of CS08, and is a precise formulation of the idea that the frequency at which the comet is perturbed by an amount 
of order J' is calculated by J'TZ{J') ~ nvb^ , where n is the number density of perturbers, v is the velocity at which 
they encounter the sun-comet system, and is the cross-sectional area for such an encounter. In words. Equation [4] 
integrates over the entire parameter space of the encounter geometry (vp,b,to, and mp), weights the integral by the 
probability density of each parameter, and uses the delta function of | A J(vp, b, io, "ip)! to select those geometries that 
produce a perturbation of size J' . 
The frequency of perturbations is linked to the distribution function through a Boltzmann equation: 



df{J,t) 
dt 



JpiJ')[f{\J' + 3\)~f{J)]d'3'. (5) 



As in CS08, the function p(J') describes the frequency per unit angular momentum space (rf^J') at which a comet 
with angular momentum J is perturbed to J -I- J'; this is the PDF of J'. We expect this frequency to depend only on 
the magnitude of the perturbation and not the direction, p(J') = p{J'), for isotropic perturbers. It is related to Tl{J') 
by integrating p( J') over the angular component of J', 'R-{J') = 2itJ'p{J'). 

We assume that the stellar perturbers have only one mass, mp, and one velocity, Vp, that can point in any direction. 
The calculation of p{J') then proceeds similarly to the calculation presented in CS08. Since the angular momentum 
excited by a perturber is proportional to nip, Vp, and b in all the same ways as the excitation of eccentricity in a nearly 
circular binary, J' cx rup /{vp9) from Equations m and H it follows that J'TZ{J') cx J'"\ and p{J') cx J'^^. 

The full calculation of p( J') requires choosing the correct expression for AJ given the timescale of the encounters. 
In the extremely non- impulsive regime (Equation [3]) , A J is averaged over ri,{t) before being used in Equation |4l For 
the impulsive case, AJ(to) retains its dependence on the position of the comet, but the subsequent integral over to in 
Equation U averages the contribution of perturbers from all possible ri,. Ultimately we arrive at the same p{J') for 
both non-impulsive and very impulsive perturbations: 

piJ')=0.74Gpa^j^, (6) 

where p = nnip, the volumetric mass density of the perturbers in space. As noted in C S08, this form of p(J') reveals 
that the angular momentum of the comet follows a Levy flight (|Shlesinger et al.lll995l) . The distribution function is 
then: 

fiJ,t) = ^^^(1 + {J/Mt)f)-'^'- (7) 

This function is self-similar, meaning that it always has the same shape centered around a characteristic angular 
momentum scale, Jc{t), that changes with time. We have chosen the normalization such that J f{J,t)d^J = 1 at all 

times. The characteristic angular momentum is near the median of the distribution, Jmodian — VSJcit)- Since the 
probability of finding the comet with an angular momentum of order J ^ Jc{t) falls off like the power law J~^, the 
mean, variance, and all higher moments of the distribution are undefined. The mean only diverges logarithmically; if 
there is a maximum angular momentum Jmax, then Jmcan = 2.3 Jc(t) log]^o(0-74 Jmax/^c(i))- 
The time derivative of Jc{t) is related to the perturbation frequency: 

Jcit) = AMGpa^ (8) 

This equation is derived by substituting the solution for f{J,t) (Equation [7]) into the Boltzmann equation (Equation 
[5]). Equation [5] determines Jc(i) even if the parameters of the perturbing swarm (p) or the comet (a) are changing 
with time. During the formation of the Oort cloud, the semimajor axes of the comets evolve as the ice giants deliver 
orbital energy to them over many interac tions. Ad ditionally, a time-varying density of perturbers may be relevant if 
the Sun formed in a dense cluster (lFernande2j|l997|). The high eccentricity but high periapse orbit of Sedna may impl y 
that the Sun was born in such an environment ([Morbidelli &: Leviso"iill2004l : iBrasser et al.|[2006t iKaib &: Quinn|[2008l) . 
A realistic statistical description of the formation of the Oort cloud must incorporate the evolution of p and a of the 
comets. 

To provide the following simple numerical example, we assume a constant p and a. The angular momentum distri- 
bution function in this case grows linearly with time, Jc{t) = 4.66Gpa^t, for Jc{t) ^ Jc{t = 0) . Using values relevant 
for the Oort cloud, we find 



Jcirc ' \0.1Mqpc-^ J KlO'^AUJ VlGyr 

where we have scaled Jc{t) by the angular momentum per unit mass of a circular orbit, Jcirc ~ \J GMqu, to make it 
dimensionless. Since our derivations neglect the non-radial motion of the comet's evolving orbit, our theory is only 
quantitatively correct for J/ Jem ^ 1- 
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This mode of growth is quahtatively different from the typical diffusive random walk. The passing stars cause a 
spectrum of perturbations that occur with frequencies inversely proportional to their size {J'TZ{J') oc J'^^)- This 
power law is such that the smallest kicks cannot accumulate fast enough to affect the distribution function. For 
example, perturbations of about the same size accumulate as a normal diffusive random walk, SJ oc •\/i/^smaU>/sinair 
In that same time, however, the comet receives, on average, a single perturbation of size SJ « J^jg oc {t/tsmaiijJsniedi- 
Thus the overall growth of the angular momentum is due to the few largest perturbations that occur over a time t. 

The distribution in angular momentum (Equation [T]) can be converted to a distribution for the comet's periapse 
distance, q, using the relation for nearly radial orbits, J — ^JiGM^q: 

f{q,t)^-^{l + q/q,{t))-^'\ (10) 

where qc{t) is the characteristic periapse associated with Jc(i). We have chosen a normalization such that / /(g, t)dq — 
1. Since Jc{t) oc t, the typical periapse distance grows as t^; the timescale for a significant change in periapse then 
depends on the comet's current q. 

These derivations of the distribution of a comet's angular momentum assumed the swarm of perturbers had a single 
individual mass and single velocity. If there are other massive perturbers with rup > Mq, such as giant molecular 
clouds. Equations [7] and [9] describe the distribution when p includes all of the perturbers: p = '^nirup^i, where and 
TTip^i are the volumetric number density and masses of the ith group of perturbers. A mass spectrum that extends 
significantly below the mass of the Sun also affects the probability distribution of the perturbations. In the generalized 
case, the slope of the perturbation spectrum sets the high J power law of the distributi on function. As long as the 
exponent of J'TZ{J') is between and —2, the angular momentum follows a Levy flight (jShlesinger et al.| [l995). For 
the precise details of deriving p{J') and f{J',t) given a general mass distribution, we refer the reader to CS08. 

4. CONNECTION TO GALACTIC TIDES 

In deriving the model presented in Section [3J we have assumed that the perturbing stars are distributed isotropically 
in Vp and uniformly in impact parameter. We expect the angular momentum distribution in that scenario to be 
axisymmetric. Field stars, which are confined to a disk with a height much less than its radial dimension, do not 
have these simplifying properties. This section uses a toy model to show how an anisotropy in the angular momentum 
distribution arises from the spatial inhomogeneity of the perturbing stars, and how this is related to the angular 
momentum distribution disc ussed in section [31 

iHeisler fc Tremaind (|1986[ ) investigated the effects of the large scale potential arising from the Galactic disk. We 
reproduce their derivation of such a torque given a simple planar model of the mass distribution. We approximate the 
disk as a stack of infinitely thin, infinitely large sheets of mass. Gauss' law shows that the sheets above and below 
both the sun and the comet produce no net acceleration on the system. The sheets that pass in between the sun and 
comet however, produce a mean torque given by: 

j = -27rGp(r6-z)(rb X z), (11) 

where p is the local volumetric mass density in perturbers, and z is the unit vector normal to the disk plane. To an 
order of magnitude, this torque is the same as our Equation [8l although it is of a completely different nature. Equation 
111! describes a smooth torque in a fixed direction, while Equation |8] is the typical value of a stochastic variable drawn 
from an axisymmetric distrib ution with zero mean. 

iHeisler fc Tremaind ()1986[ ) also performed numerical experiments to verify that on very long timescales, stellar 
scattering indeed produces a mean growth on top of the stochastic e volution. The importance of the Galactic tides 
has been appreciated in su bsequent studies of Oort cloud dynamics ([Duncan et al.l [1987t [Heisleij [19901 : [Pones et al.l 
[20041 : [Rickman et al.l [2008[ ) , although the relationship between the stellar encounters and the tidal torques is rarely 
addressed. Tidal torques are usually treated as separate from the effects of stellar encounters, even though the torque 
is provided by the same stars that cause the stochastic evolution. By adapting our formalism to reflect a planar 
distribution of perturbers, we reproduce the effects of the Galactic tides, and in doing so find the distribution function 
that accounts for both modes of angular momentum growth^ 

We follow the example of the numerical experiments of [Heisler fc Tremaind ([19861 ) and approximate the Galaxy 
locally as a uniform disk of material, with a height much smaller than the scale of the other two dimensions. To 
create the planar symmetry in the model of stellar encounters, the velocities of the perturbers are restricted to a 
single direction. While this is not a realistic representation of the directional distribution of field star velocities, it is 
a simple model to explore and provides a clear example with which to examine the effects of a velocity asymmetry. 
With Vp fixed, the impact parameter b is confined to a plane, the aspect ratio of which has a much smaller height 
than width. Both of these properties, a single direction for Vp and a non-unity aspect ratio, introduce asymmetries in 
the distribution function of the comet's angular momentum. 

For isotropic perturbers, perturbations of any size J' occur with the same likelihood in all directions in the plane 
perpendicular to rf,. This ensures that the mean of 3{t) is zero, even though the typical magnitude of the angular 
momentum increases linearly with time. The cross-section for an interaction in the tidal limit {b ^ r^) scales as 5^, 
which fixes the power law of the single perturbation PDF. In the planar model, the cross-sectional area that contributes 
perturbations with small J' is less than for impact parameters larger than the disk height. The contributions of these 
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Fig. 1. — Contours of constant J' on the space of impact parameters b/(ri, ■ z) for positive and negative values of each component of the 
vector perturbation. The levels are spaced in multiples of ten from J' /{2Gmp/vp) = ±10~* to ±1. Panel b, which shows the contours for 
the y direction, is symmetric with respect to positive and negative perturbations. The center of panel a shows an isolated region of negative 
X perturbations that causes an asymmetry in the distribution function. By randomly sampling this space of impact parameters we generate 
the PDF of the perturbations. The marginal PDF for positive and negative J'^ are plotted in panel c; the spike contains perturbations 
from the central region of panel a and is the source of the Galactic tidal torques on the comet. 

regions to each component of J' depends on the angle between the comet and the disk plane so the axisymmetry is 
broken. However, these differences manifest only in the lowest J', and their effects on the distribution of accumulated 
angular momentum are always washed out by the larger perturbations from impact parameters less than the disk 
height. 

Another asymmetry results from the impact parameters of 6 ^ r^. For b > rb, there is as much cross-sectional area 
contributing positively to each component as there is negatively. Impact parameters that pass between the sun and the 
comet, however, impart angular momentum in one direction of one component only, depending on the angle between 
ff, and Vp. Not coincidentally, the mean torque found in the smooth distribution limit, Equation is attributed to 
the disk of stars passing between the Sun and the comet. 

We quantify the effect of this asymmetry by calculating the marginal probability density of each component of 
the angular momentum vector due to single interactions. Since we have lost the symmetry that admitted the simple 
analytic solutions, we employ a Monte-Carlo procedure. The position of the comet, which we hold fixed in this example, 
is Tf, = y + z, so the Sun-comet distance is rb = a/2. The perturber velocities are set to the z direction: Vp — — z. 
The possible impact parameters of the perturbers are then restricted to the x — y plane. We randomly choose impact 
parameters such that they are uniformly distributed over the plane and calculate the A J delivered to the comet. We 
assume the other parameters of the system are held constant {vp and m^), and to reduce the notation, we use units 
where IGmpjvp = 1. The angular momentum is confined to the plane perpendicular to rb, which in these coordinates 
is defined by the basis vectors x and {y — z) / \/2. For simplicity we discuss the x and y components of the perturbation, 
A J ■ X — J'^ and A J - y — J'y. In the z-direction, A J • z is exactly the same as J'y. The positive and negative values for 
J'^ and J' are binned separately; the resulting four histograms then describe the marginal PDF for each component. 

Figure [T] illustrates the calculation of the single interaction PDF. Panels a and b show logarithmically spaced 
contours of constant J'^ and J'y respectively in the plane of possible impact parameters, with the other parameters 
of the interaction fixed (r;,, Vp, m;,). The impact parameter plotted is scaled by • z = = 1. The solid contours 
correspond to positive perturbations and the dashed lines to negative ones. In panel b), the contours for ± exhibit 
an axisymmetric pattern; for each unit of area that contributes perturbations of a given magnitude greater than 
zero, there is an equivalent area where perturbations have the opposite sign. Thus the single interaction marginal 
PDF of perturbations in the y directions are identical and unchanged from the isotropic case: J'y ^^ for the distant 
perturbations, J'y{b S> ^b), and J'y~^ for the close encounters, Jy{b ^ Tb). There is no coherent accumulation of 
angular momentum in the y direction. 

The contours of panel a), while symmetric at larger b, are not symmetric in the center, where the perturbations 
only add angular momentum in the negative x direction. There is no equivalent area that delivers angular momentum 
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Fig. 2. — Marginal distribution functions of two components of the angular momentum as a function of time. The dotted and dot-dashed 
lines plot the marginal distribution, dN/d{log \Jy\), of the y component of the angular momentum J{t). The thick line is the distribution 
of the X component when it is negative, and the dashed line is the positive side. In the top two panels, the comet's angular momentum is 
best described by the Levy flight behavior caused by stochastic stellar perturbations. In the bottom two, the coherent torque attributed 
to the Galactic tides dominates the evolution, causing a visibly asymmetric distribution. The thin line in the bottom panel is a Gaussian 
distribution with the mean given by the Galactic tidal torques and the variance given by the variance of the single interaction PDF 
multiplied by the number of encounters. 

with the opposite sign. We plot the marginal PDF of J^, \ J^\TZ{J!^), in panel c) of Figure[l] where the solid line is for 
perturbations where > and the dashed line is for < 0. The values along the ordinate represent the probability 
of perturbations with strength of order relative to the lowest value plotted. In the tidal and close encounter regimes, 
the two functions are identical. For of order unity, the contribution of the central region in panel a) is obvious. It 
is these interactions that give rise to the torque associated with the Galactic tides. 

The marginal PDF of highlights the source of the Galactic tidal torque. However, it remains to describe how this 
manifests in the time-dependent distribution function of the comet's angular momentum. In Section |31 we used the 
Boltzmann equation (Equation [5]) to relate the axisymmetric single perturbation PDF (p(J')) to the distribution of 
angular momentum (/(J(t))). That derivation, however, depends on the simplifications afforded by the single power 
law form of p{J')- For the non-axisymmetric single perturbation PDF depicted in Figure [l] an analytic solution to 
the corresponding Boltzmann equation would be much more difficult to calculate. 

Instead, we use a bootstrap technique to estimate the distribution function from a sample of single perturbations. 
The velocity of the perturbers, Wp, their number density, n, and the area sampled when generating the single interaction 
PDF, TT&^j^jj, set the average time associated with each perturbation, l/r = mrb'^^^Vp. The angular momentum at a 
time t is then the sum of t/r single perturbations. By randomly choosing t/r perturbations from the PDF and adding 
them vectorially, we generate a sample of angular momentum vectors that reflect the distribution function at that 
time t. 

To accurately probe the evolution over many orders of magnitude, several single interaction PDFs with different 
^max were used. Ignoring large impact parameters increases t, or equivalently, samples the close encounters more often 
over a fixed number of perturbations. We verified that the distribution functions calculated with large r (small 6max) 
are not significantly affected by ignoring the frequent perturbations of smaller J'. 

The marginal distribution functions at four different times are shown in Figure [2] Each histogram contains 10^ 
bootstrapped J{t), generated from the sum of between 4 and 1000 single perturbations. The distribution of Jy{t) is 
plotted in the dotted lines for Jy{t) > and dash-dotted for Jy{t) < 0. For Jx{t), the solid line represents the negative 
perturbations and the dashed line the positive ones. 

The top panel shows the angular momentum distribution at early times, or equivalently, at low typical angular 
momenta. For reference, we denote this time ^o- Since the single interaction PDF for perturbations of this magnitude 
is axisymmetric, all four functions are identical. The excess of perturbations to negative is not visible as the 
likelihood for those encounters is too low to be sampled in the 10^ vectors generated for the plot. 

The second panel depicts the four distribution functions 100 times later than the time of the top panel. Again 
both functions show a similar shape, and the typical value for all four has grown linearly with time as predicted by 
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Equation [8] The trajectories passing between the sun and the comet have been sampled in a smaU fraction of the 
generated 3{t), and the contribution from the spike of Figure [1]; is apparent. AdditionaUy the normahzation of the 
positive distribution of Jxit) has fallen to reflect the breaking of the symmetry around = 0. The distributions 
in the first and second panel can be said to be dominated by the influence of the stellar perturbations, and are not 
strongly affected by Galactic tides. Although the mean of the distribution is always set by the tides (see Equation 
here this value of angular momentum is only realized after rare but strong interactions. The most likely angular 
momentum vectors, at early times, are distributed axisymmetrically around the origin. 

In the third panel the non-axisymmetric growth is manifest. Due to the higher slope of the single encounter PDF, 
the distribution of the y component of the angular momentum has begun to grow only as i^/^; the accumulations of 
kicks from all of the impact parameters smaller than ri, contribute to the shape of this distribution. Unfortunately 
a PDF of this slope does not admit a self-similar distribut ion function; asymptotically, the distribution approaches a 
Gaussian logarithmically over time ([Shlesinger et al.lll995l) . 

The perturbers passing between the sun and the comet deliver angular momentum in the —x direction coherently 
and thus the typical —Jx{t) continues to increase linearly in time. The normalization of the histogram for positive 
Jx{t) has decreased substantially, which is another indicator that the total distribution of Jx{t) is no longer centered 
on the origin. In the fourth panel, only 10 times later than the third, the marginal distribution function for Jx{t) 
is entirely dominated by the accumulated effects of non-canceled encounters. There are no values of Jx{t) > in 
the sample at this time. Again, the distribution function does not admit an analytic form. For reference, we plot a 
Gaussian distribution with the mean described by Equation 1111 and the variance expected given the single encounter 
PDF, = (TpDpi. The distribution function only approaches this approximated shape logarithmically in time. 

Figure [2] reveals the nature of the coherent torque by Galactic tides as merely the long term effects of anisotropic 
stellar encounters. It is only a matter of principle what to call the interactions of the comets with field stars. To 
determine the relevant behavior, one must specify which impact parameters are the most important for the behavior of 
the comet. On shorter timescales, or for smaller angular momenta, the distant perturbations create the axisymmetric 
distribution function associated with stochastic stellar encounters. Over timescales long enough that many trajectories 
have sampled the region between the Sun and comet, the system is best characterized as evolving under the Galactic 
tides. 

As a physical example, we again examine the formation of the Oort cloud, where a proto-comet must gain enough 
angular momentum to raise its periapse q by Ag to avoid perturbations from the planets. The influence of the p lanets 
falls off rapidly with increasing q, so a reasonable value for Aq/q is on the order of 10% (jO uncan "eTan fToSTt). The 
distant stellar encounters will be responsible for building the Oort cloud if a single interaction at an impact parameter 
6 ~ a can provide enough angular momentum to increase the periapse. If these single encounters are too weak, the 
coherent growth due to Galactic tides is required. We find the following inequality for when the mean tidal growth, 
rather than stochastic evolution, dominates: 




where Vq = (Grrio/g)^/^ is the local rotational velocity at periapse. At the semimajor axis of Jupiter, this velocity is 
about 15 km s^^, and near Neptune it is about 5 km s~^. Typical velocity dispersions of stars in the solar neighborhood 
are 15 — 40 km s^^ (Binney & Trcmainc 1987). Then in the inner solar system, the tidal torque is less important than 
the stellar encounters for freeing the comets from planetary perturbations. In the outer solar system, the left hand 
side of Equation [1^ is close to unity, meaning the stellar encounters and the tidal torque play a comparable role. 

Our new understanding of the relationship between stellar encounters and tides presents a clearer picture of the most 
appropriate way to model the excitation of angular momentum in an Oort cloud comet. If the prescription for stellar 
encounters includes the planar symmetry of the stars, then no extra torque is required to represent the Galactic tides. 
If the stellar encounter model has an isotropic velocity distribution, then an extra term representing the torque should 
be included, but only at late enough times that encounters passing between the sun and the comet are common. 

5. CONCLUSIONS 

In this work we have shown that the angular momentum delivered to nearly radial comets by passing stars follows 
a Levy flight. From the properties of a single scattering between the comet and the star, we derive the distribution 
function of the angular momentum of the comet as a function of time. Our calculations agree with the estimates made 
in earlier work on Oort cloud formation, that stellar perturbations can raise the periapses of comets significantly in only 
several hundred Myrs. A careful examination of the scattering process for an anisotropic velocity distribution reveals 
the presence of the coherent angular momentum growth that is usually attributed to the large scale potential of the 
Galaxy. The effects of stellar encounters and the Galactic tidal torques then cannot be treated as two distinct processes. 
On shorter timescales the distribution function of the comet is unaffected by the tidal torque; on long timescales the 
distribution is entirely dominated by it. Since the presence of the tidal torque depends on the perturber velocity 
distribution, simulations of cometary evolution that include stellar encounters must be careful not to double-count the 
Galactic tides by either including an explicit torque or enforcing a planar symmetry, but not both. 

These results provide a formal understanding of the effects of stellar encounters on nearly radial comets, but it is only 
the first step towards a complete statistical picture of the formation of the Oort cloud. The shape of the distribution 
function of the angular momentum at early times will not be entirely isotropic due to the triaxial velocity distribution 
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of field stars; however, tliis anisotropy will be overwhelmed at the current epoch by effects of the Galactic tidal torque. 
The effects of the stellar perturbations must be convolved with the diffusion of the comets' semimajor axes caused 
by planetary perturbations. This type of diffusion is not without complications, as o rbital resonances bet ween the 
comet and the pl anet must be accounted for to produce accurate diffusion coefficients ([Malvshkin fc Trem ainc 199§[, 
iPan fc Sarill2004l ). Additionally, the diffusion of the semi major axis for a comet whose orbit crosses that of a planet 
has been shown to exhibit properties of a Levy flight (Zhou et al.ll2002^ . 

We thank the Institute for Advanced Study for their hospitality while some of this work was completed. R.S. is a 
Packard Fellow. This work was partially supported by the ERG. 
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